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ABSTRACT 

Working in a spherically symmetric ansatz in Minkowski space we discover new solutions 
to the classical equations of motion of pure SU (2) gauge theory. These solutions represent 
spherical shells of energy which at early times move inward near the speed of light, excite the 
region of space around the origin at intermediate times and move outward at late times. The 
solutions change the winding number in bounded regions centered at the origin by non-integer 
amounts. They also produce non-integer topological charge in these regions. We show that 
the previously discovered solutions of Liischer and Schechter also have these properties. 
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I. INTRODUCTION 

In non-Abelian gauge theories coupled to fermions, gauge field configurations with non- 
trivial topological charge can cause explicit violation of conservation laws as a result of anoma- 
lous Ward identities.-^ This leads to the violation of chirality in QCD and non-conservation 
of baryon and lepton number in the electroweak theory.^ The traditional approach to un- 
derstanding these effects uses semiclassical barrier penetration where the tunneling solutions 
are Euclidean instantons.^ In the electroweak theory the height of the barrier is of order 
M^ja^ ~ lOTeV - the energy of the sphaleron^ - and at energies much below this, fermion 
number violation is exponentially suppressed. 

More recently it has been suggested that fermion number violating processes may become 
unsuppressed if high temperatures^ or high energies® are involved. The idea is that instead 
of tunneling through the barrier the field configuration can pass over it. 

At temperatures above the barrier height it is generally assumed that field configurations 
which pass over the barrier can be found in the thermal ensemble. The rate estimates proceed 
through thermodynamic arguments which are fairly general and make reference to little more 
than the temperature, barrier height and coupling constant. 

At high energies the existing techniques are based mostly on Euclidean methods which 
deal with contributions of instanton-like configurations to the functional integral (for recent 
reviews see Ref. [7]). These essentially semiclassical methods have intrinsic problems at en- 
ergies of order the barrier height, exactly where anomalous fermion number production may 
become unsuppressed. 

We are pursuing a complementary approach^ based entirely in Minkowski space which 
corresponds to passage over the barrier. Our approach may be separated into two main parts: 
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the creation and evolution of finite energy classical gauge fields by particle collisions and the 
dynamics of fermion number (or chirality) violation in the presence of such background gauge 
sources. In this paper we study the classical evolution of certain gauge field configurations in 
Minkowski space. The problem of fermion production is addressed in a companion paper. ^ 

In the Euclidean approach what is typically done is to only consider finite action contribu- 
tions to the functional integral. These configurations fall into homotopy classes distinguished 
by an integer, Q, called the topological charge, given by 



What we discuss here are Minkowski space classical solutions that produce arbitrary non- 
integer topological charge and winding number change. 

Consider a spherical region of space with radius R containing vanishing gauge fields at 
some large negative time, —Tq. Imagine sending in a localized spherical wave- front of energy 
E. At time t ~ —R the front first reaches the region. Energy flows in, bounces back at time 
t ~ 0, and then leaves the region at time t ~ i?. At a large positive time, -I-Tq, the region 
of space interior to R is once again in a pure gauge, but now the vector potential need 
not vanish. The question we address is, what is the topological structure of the pure gauge 
configuration left behind by the wave- front? We will construct solutions to the Yang-Mills 
equations which describe such localized wave-fronts and study in detail the winding number 
of the pure gauge configurations that the fronts leave behind. In addition we compute the 
topological charge created inside the region. We will also study the topological properties of 
the previously discovered solutions of Liischer^^ and Schechter,^^ which we also show describe 
localized incoming and outgoing spherical waves. 




(1.1) 
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Classical vacuum configurations are typically classified-*^ by the winding number of the 
associated gauge function. In the = gauge the classical vacuum configurations are of the 
form Ai — i/g UdiU^ where t^(x), the gauge function, is a map from the into the gauge 
group. The winding number is 



gauge group. These maps can be classified by an integer which labels the homotopy class of 
the map. Formula (1.2) is an explicit expression for this integer. The standard picture^^ of 
the gauge theory vacuum is based upon this classification. The ^-vacua of QCD are linear 
superpositions of states associated with classical vacuum configurations which have integer 
winding number. However, if ?7(x) is not required to approach a constant at spatial infinity, 
such an elegant vacuum classification is not available. One typically argues, however, that 
by causality the values of fields at spatial infinity are unimportant, so imposing a convenient 
boundary condition does not entail a loss of generality. Also, instantons can be viewed as 
changing winding number by an integer amount. If only instanton-like configurations are 
important, then restricting to vacua with integer winding number can be justified after the 
fact. 

However, if we return to the description of the spherical region of space with radius R 
that is excited by a wave-front and then relaxes back to pure gauge, we find that U ai R does 
not equal a direction-independent constant. Accordingly, the local winding number at late 
times, computed by doing the integral in (1.2) only out to R, does not equal an integer even 
if R is very large. That is to say, if our known universe is the region in question, "aliens" 





If we assume that U{x) 



1 as I x 1^ oo then U{'x.) is actually a map from into the 
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from the outside could send in a wave which could leave our universe in a different classical 
vacuum, one which did not belong to the previous topological classification. 

If the wave which visits the region changes its winding number by a non- integer amount, 
and creates non-integer topological charge, we can ask what the implications are for anomalous 
fermion production in the region. In a companion paper^ we study fermion number production 
in a (l+l)-dimensional theory with an arbitrary background gauge field. We arrange for the 
background field to change its winding number in local regions and we allow for these local 
changes to be non-integer. We find that fermions are produced only in these local regions. 
Furthermore we discover that the the number of fermions produced in a local region, in a 
quantum average sense, is equal to the change in winding number of the gauge field in this 
region. That is, if in a local region the change in winding number is /, and the background 
field is used time and time again to quantum-mechanically produce fermions, then in the local 
region with each observation you will find integer net fermion production, but the quantum 
average will be /. This suggests that if we restrict our attention to the region of radius 
R, and a spherical non-Abelian wave comes in and leaves, there will be fractional fermion 
number production in the sense described above. Here we are discussing net fermion number 
production in the region, that is the difference between the number of fermions coming into 
the region through the two-sphere of radius R and the number leaving. 

The outline of the paper is as follows. In Section II we explain the spherical ansatz"^^ which 
we use to simplify the Yang-Mills equations and we introduce a further ansatz which allows us 
to find new solutions. In Section III we study in our language the solutions previously found 
by Luscher^° and Schechter.^^ We show that these solutions have many of the qualitative 
features of the solutions of Section II. In both cases the solutions have certain transformation 
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properties under S0{2, 1) which allows us to generate new solutions. The action of symmetry 
transformations on the solutions is discussed in Section IV. In Section V we investigate the 
topological properties of both classes of solutions. Finally in Section VI we discuss further 
questions which arise from this work. 

II. THE SPHERICAL ANSATZ AND NEW YANG-MILLS SOLUTIONS 

In this section we consider the spherical ansatz^^ for SU (2) gauge theory. We follow the 
notation of Ref. [14]. We find new classical solutions in (3+l)-dimensional Minkowski space 
and discuss their properties. 

The action for pure SU (2) gauge theory is 



where F^i, = F^^{a°'/2) = d^A^, — dj^A^ — ig[A^^Aj,] is the field strength and A^ = 
A'^^{a°'/2). Throughout we use the space- like convention for the metric tensor, ?7^j^ = 
diag (—1, +1, +1, +1), The spherical ansatz is given in terms of the four functions ao, ai, a, (3 




(2.1) 



by 



^o(x,t) = — ao(r,t)o--x 




(2.2) 



where the matrix- valued functions {e^} are defined as 





i [a ■ xcTj - Xi] = eijkXjak , 



(2.3) 



4 = 



and where x is a unit three-vector in the radial direction. 
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The action (2.1) in the spherical ansatz takes the form 

^ = ^/^^/°°^-(-^-VM.r^-(i^Mxr^^x-^(ixi^-i)') • (2.4) 

The (1 + 1) -dimensional field strength is defined as f^i, = d^ai, — diya^, where jJi-iV = t,r, 
and indices are raised and lowered with fj^jy = diag(— 1, +1). A scalar field x = cu + with 
covariant derivative D^x = {d^ — ia^) x, has also been introduced. 

The ansatz (2.2) preserves a residual U{1) subgroup of the SU{2) gauge group consisting 
of the transformations, 



U (x, t) = exp 



These induce the gauge transformations 



in{r, t) 



a ■ X 



(2.5) 



O/i ^ + ' X exp(il])x , 



(2.6) 



which are seen to leave (2.4) invariant. 

Vacuum configurations in (3 + l)-dimensions in the = gauge are given by the pure 



gauges, 



(2.7) 



where C/ is a time- independent function from three-space into SU (2) . If we wish to maintain 
spherical symmetry then 



U (r) = exp 



a ■ X 



(2.8) 



In terms of the reduced theory we have from (2.7) and (2.2) that 



dn , 

J ' X 
ar 



= —le 



(2.9) 



This form for the vacuum configurations in the (l+l)-dimensional theory could easily be 
obtained from the action (2.4). 

In (3+1) dimensions the function (2.8) is well-defined at r = only if f2(0) = 27rn with 
n an integer. If the further restriction is imposed that U ^ 1 ai spatial infinity, then we have 
n(oo) = 2nm with m an integer. The winding number of this configuration, as computed by 
(1.2), is m — n. If we impose no boundary condition at spatial infinity, then U would not 
have an integer winding number. From the point of view of the reduced theory, the vacuum 
configuration is given by (2.9) and we see no reason, at this point, to impose restrictions on 
Q at zero or infinity. We will see, however, that a boundary condition at r = (but not at 
r — oo) will be dynamically imposed on all finite energy solutions to the equations of motion. 

The (1 + 1) -dimensional equations of motion for the reduced theory (2.4) are given by 

-9^ (rVp.) = i [{D^xT X - X*D,X] , (2.10a) 

(^-D2^^(|x|'-1))x = 0. (2.10b) 
Let us express the complex scalar field x polar form, 

X{r, t) = -ip{r, t) exp [i(fi{r, t)] , (2.11) 

where p and (p are real scalar fields. We must be careful in (2.11) when p vanishes, since the 
angle (p can change discontinuously by tt if we require that p > 0. We will choose the sign of 
p to ensure the continuity of <^ as x varies smoothly through zero. In terms of p, (f and a^, 
the four equations contained in (2.10) read 
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d''{r^f^,)+2p^d,<p-a,) = , (2.12a) 

di'd^p - p {d'^v - a'') {d^<fi - a^) - (p2 - 1) = , (2.12b) 

and 

d^'lp^d^ip-a^)] = . (2.12c) 

The last equation follows from (2.12a) so we see there are only three independent equations, 
not four, as we expect because of gauge invariance. Since in (1+1) dimensions /^j^ must be 
proportional to e^y, it is convenient to define a new field ip by 

r^U, = -2e^,i^ , (2.13) 

where e^j^ is the antisymmetric symbol with eoi — +1. Contracting (2.13) with €^^'^ gives the 
relation 

^'^'^V-^^V' , (2.14) 
which we use below. By using (2.13) in (2.12a) and contracting with e^"' we find 

d<^^ = -e^- p\d,ip - a,) . (2.15) 

Notice that (2.15) implies (2.12c). Prom this last equation it follows that 

^4^) = ^"'^«a.= ^V', (2.16) 

where we have used (2.14) in passing to the last equality. This gives an equation solely in 
terms of the fields p and ip. We may also use (2.15) to express the second term in (2.12b) in 
terms of only p and ip. We then have the alternate classical equations of motion 
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-djp + dip - 1 {dt^f + 4 - ^P^P" - 1) = , (2.17a) 

pC j.^ 



Equations (2.17) are equivalent to the original Eqs. (2.10), but for our purposes they will 
be more convenient. Note that p and if: are gauge invariant fields, and we have only two 
equations. 

Using the equations of motion, the energy associated with the action (2.4) can be written 
in terms of p and i/^ as 
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We are interested in finding finite energy solutions to (2.17). Thus the form of p and ip when 
r ~ must be 

p = 1 + 0{r^/^ + ^) , (2.19a) 
V; = e»(rV2 + ^) , (2.19b) 

where A, 7 > 0. We have used the symmetry p — > —p under which (2.17) is invariant to fix p 
to be at +1 and not — 1 as r goes to zero. 

Witten^"^ observed that (2.4) is the action for an Abelian Higgs model on a surface of 
constant curvature. To see this in terms of p and fields we can write the equations of motion 
(2.17) in covariant form 
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^ 9m iV9 g^^d.p) + 4 Q^^'d^^d-^ - PiP^ - 1) = . (2.20a) 
y9 P 

^ 5^ [yg 9^"^^ - 2V; = , (2.20b) 

where y^i^ = fj^^/r'^ and (7 = |det(7^,^|. (Recall that the covariant derivative of a vector 
field can be written as V^V^ = g~^^'^ dfj,{g^/'^ V'^).) The metric g^i^ is the metric for two 
dimensional De Sitter space. To see this, consider the hyperboloid Zq — zj — Z2 = —1 where 
the Zi are parametrized in terms of r and t as 

zo = , 2.21a 

zr 

zi = - , (2,21b) 
r 

It is straightforward to check that ds'^ = dzQ — dzl — dz^ = {—dt'^ + dr'^)Jr'^. The coordinates 
r and t cover only half of the hyperboloid, , for which Z0 + Z2 > 0. This is shown in Fig. 1. 
The ^-coordinates (2.21) will be useful for studying symmetry properties of solutions. To 
obtain our solutions, however, it is more convenient to work with coordinates w and r that 
live on the hyperboloid itself, and in fact cover the whole hyperboloid: 

zq = — tanw , (2.22a) 
zi = sinr/ costu , (2.22b) 
Z2 = cost/ cos , (2.22c) 

where |t| < tt and \w\ < 7r/2. The coordinate is a bounded measure of the longitudinal 
position along the hyperboloid and r measures the azimuthal angle. Notice that for fixed t, 
as r varies from zero to infinity, w varies from — 7r/2 to +tt/2. In terms of w-t variables, the 
metric is ds^ = {—dr"^ + dw'^) / cos^ w, so (2.20) takes the form 
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-dip + alp-^ (a^i, f + 1 (S„^)^ - = , (2.23a) 

-Bj?4)+aJ?4)-^ = 0. (2.23b) 

\ P J \ P J cos^w 

We make the ansatz that there are finite energy solutions to (2.23) which are independent 
of r. These obey the ordinary difi'erential equations, 

^,^(^y_^!0^ = 0, (2.24a) 

aw^ p-^ \ aw J cos^ w 

d f 1 dilj\ 2ip 



A \ 2 ^ t 2=0- (2-24b) 

aw \ aw J cos^ w 

It is interesting to note that the previously discovered solution of de Alfaro, Fubini and 

Furlan^^ can immediately be obtained in the following form: 

^ = , p= -sinw . (2.25) 

Before we consider the full coupled system (2.24), we examine the ■0 = sector where 
there are solutions other than (2.25). Equation (2.24a) with '0 = describes a "particle" 
rolling in a "time" dependent potential V = — (p^ — 1)^/4 cos^ tu. We can make a change of 
variables from w to rj, defined by sinw = tanh?/, which removes the "time" dependence from 
the potential at the price of introducing a "friction" term: 

p + tanh?7p - p(p^ -1) = . (2.26) 

The dot denotes differentiation with respect to r], and r] runs from — oo to +oo. Now (2.26) has 
a mechanical analogue to a "particle" rolling in a potential U = — ;|(p^ — 1)^ with a "frictional" 
force — tanh?7 p. The "energy" of the mechanical system is defined by T = ^p'^ + U, and for a 
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particle satisfying the equations of motions, T = — tanhr/ p^. So for negative r] the frictional 
force pumps energy into the system and for positive rj it pumps energy out of the system. We 
are interested in starting the particle at p(— oo) = 1 because of (2.19a). The particle can only 
stop rolling at an extremum of the potential U so that p(+oo) = ±1,0. If the particle does 
not stop rolling, p grows without bound which is unacceptable if we wish (2.18) to remain 
finite. If the energy pumped into the system at negative 77 exactly matches the dissipation 
at positive rj, then p can climb up the other hump at p = —1, as does the solution (2.25). 
Typically this does not happen. If the velocity of the particle when it reaches the bottom of 
the potential is too large the particle overshoots the other hump and p blows up at finite 77. 
However, if the velocity of the particle is not too large, it gets trapped at the bottom of the 
potential and dissipation drives p to zero. An example of a numerically generated solution 
to (2.26) which goes from p = 1 to p = is shown in Fig. 2. We can specify a general 
trajectory by choosing p and p at ry = 0, giving a two parameter family of solutions. For a 
given p(0), we must fine-tune p(0) so that p{r]) ^ 1 as 77 — > —00. This freezes out a degree 
of freedom and gives a one parameter family of solutions. The future behavior is typically to 
have p dynamically driven to zero. However with further fine-tuning we can arrange for p to 
approach —1 at late times, which is the solution (2.25). 

At first sight making the r-independent ansatz for (2.23) looks similar to making the 
assumption that solutions to (2.17) exist which are independent of time. It is known that 
such static solutions do not exist. We illustrate this by considering (2.17a) with '0 = and 
assuming dp/dt = 0. Making the change of variables s = Inr, the equation of motion for p 
becomes 

p" - p' - p{p^ - 1) = , (2.27) 
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where the prime denotes differentiation with respect to s. This represents a particle roUing 
in an inverted potential U — —\{p^ ~ 1)^ with a constant frictional force +p'. For a particle 
satisfying the equations of motion, the power input is T' = p'^. Thus energy is always pumped 
into the system, and a particle that starts at p = 1 will always overshoot the other hump. 
Thus, even though r and t are both time-like, making the static r ansatz in (2.24) is in fact 
quite different from making the static t ansatz in (2.17). 

We are now interested in finding finite energy solutions to (2.24) with non-trivial ip- 
dependence. Recall that — 7r/2 < w < 7r/2 and that for fixed t as r — > 0, — > — 7r/2 and 
as r — > oo we have w — > 7r/2. In the appendix we investigate the form of finite energy 
solutions near the end points of the range of w. We find the following asymptotic behavior 
for w = -7r/2 -|- e with < e « 1, 

p=l + 0{e^) = l + 0{r^) , (2.28a) 
ijj = 0{e^) = 0{r^) . (2.286) 

Note that these conditions, imposed by dynamics, are stronger than the finite energy condi- 
tions (2.19) . 

We also show that for finite energy solutions in which ip is not identically zero, p is 
dynamically driven to zero and undergoes an infinite number of oscillations as w — > n/2. 
Writing w = tt /2 — e with < e << 1, the form of such solutions to leading order in e is 

p = Ae^/^ cos[clne + 5] , (2.29) 

where c = -\/3~+T6^/2 with ipQ = 'ip{T:/2) and A, B are integration constants. 

We numerically found solutions that are consistent with the above asymptotic limits 
using a Runge-Kutta method. We integrated equations (2.24) starting at a point slightly to 
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the right of — 7r/2 with the initial condition that p = +1 and = 0. We were then free to vary 
the initial derivatives of p{w) and '^(it;) giving a two parameter family of solutions. Typical 
solutions are shown in Figs. 3-4. Note that p approaches zero in an oscillatory manner at 
the right end point, consistent with the above discussion. 

We can take the numerically generated solutions to (2.24) in terms of w, and using (2.21a) 
and (2.22a) convert them into space-time solutions with r-t dependence. For very early and 
late times these solutions propagate undistorted keeping their shape. We may substitute 
these solutions into the integrand of (2.18) and obtain an energy density profile. Fig. 5 shows 
this profile for a sequence of times. The energy density forms a localized shell which moves 
undistorted in a soliton-like manner at early and late times. This early and late time behavior 
can be understood since our solutions depend only on zq given by (2.10a). For example, for 
t » 1 and r near t, zq approaches t — r, which shows that our solutions move undistorted 
near the speed of light. The behavior of the energy shell as a function of time is consistent 
with the observation of Coleman and Smarr that the radius of gyration, f(t), for a shell of 
pure glue must obey f(t)^ = t"^ + where tq is a constant. 

We will see in Section V that only the solutions with ip ^ have interesting topological 
properties. However, before studying the topological structure of these solutions, we exhibit 
in the next section a class of previously discovered solutions whose topological properties we 
will also investigate. 

III. THE SOLUTIONS OF LUSCHER AND SCHECHTER 

There are other classical solutions in pure SU (2) gauge theory discovered by Luscher^° 
and Schechter.^^ We will show that these solutions can also be described as moving spherical 
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shells of energy. Motivated by the conformal invariance of Euclidean solutions, Luscher and 
Schechter constructed an SO {4) symmetric Minkowski solution, where the SO {4) is a subgroup 
of the (3+l)-dimensional conformal group <S'0(4, 2). This 50(4) contains the rotation group, 
so their solution can be cast in the spherical ansatz. Converting their solution to the notation 
of the previous section we find that a^, a and /3 are expressed in terms of a single function 
^(t) as 

Op = -q{^) d^w , 

a = q{r) sintu costu , (3-1) 

(3= -{1 + q{T) cos^ w) , 
where /j, — t,r. Expressing p and •0, which now depend on both r and w, in terms of (3.1) we 

find 

P^{w,t) = l + q{q + 2)cos^w , (3.2a) 
ip{w^T) = - qcos^w , (3.2b) 

where the dot denotes differentiation with respect to r. By substituting (3.1) into (2.10) or 
equivalently and more simply substituting (3.2) into (3.23), it can be seen that qir) obeys: 

g + 2Q(Q + l)(Q + 2) = . (3.3) 

Note that (3.3) is the equation for an an-harmonic oscillator with the potential XJ = 
lq^{q+2y. The solutions of (3.3) may be characterized by the "energy" e of the associ- 
ated mechanical problem, 

e^^q' + U{q) . (3.4) 
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There are two classes of solutions depending on whether e is smaller or larger than 1/2, the 
barrier height of U{q) at the unstable point g = — 1: 

q{T) = -1 ± (1 + v^)V2dn[(l + V2i)^/\T - To); fci] 

(3.5) 

kl = 2V2e/il + V2e); e<l/2 , 

and 

g(r) = -1 + (1 + v^)V2cn[(8£)i/4(T - To); ^2] 

(3.6) 

A;^ = (l + v^)/(2v^); e > 1/2 , 
where dn and cn are the Jacobi elliptic functions and To is an arbitrary parameter. For 

£ < 1/2, the "particle" can be trapped in either well, which corresponds to the two forms for 

g(T) in (3.5), while there is only one solution when e > 1/2. The parameter to corresponds 

to the time at which the particle moving in the potential U (q) with energy e is at a turning 

point. 

Note that the solutions (3.5) and (3.6) represent bound motion in a potential and are 
therefore periodic in the variable t. The period depends upon e and in general is not 2n 
divided by an integer. Therefore g(7r) does not match g(— tt) and, the fields p{w,t) and 
i(^{w,t) are generally discontinuous along the line on the hyperboloid t = ±7r. This line, 
however, lies outside the physical region . Only for special discrete values of e for which 
the period is 27r divided by an integer will p and i/j match along the t = ±7r line and therefore 
be continuous on the whole hyperboloid. 

As was shown in Ref. [10], the Luscher-Schechter solutions have finite Minkowski action. 
The solutions of Section II, however, have infinite action. 

Like the solutions of Section II, the solutions described here also give spherically symmet- 
ric waves of localized energy density. It can be shown that the energy density is independent 
of To and depends only on the parameter e. Fig. 6 shows some r-profiles of the energy density 
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for a sequence of times. Again in the distant past the solution propagates undistorted in a 
sohton-hke manner with the energy density locahzed in a spherical shell which is collapsing 
at near the speed of light. The shell becomes small, distorts, and bounces back producing an 
expanding shell. As the shell expands it leaves the region of space behind it in a pure gauge 
configuration. 

IV. SYMMETRY TRANSFORMATIONS OF THE SOLUTIONS 

In this section we show how to construct families of solutions by acting with symmetry 
transformations on the solutions we have described in the previous two sections. The reader 
who is primarily interested in the topological properties of the solutions can skip to the next 
section. 

The fields p and '0 which satisfy (2.20) live on the hyperboloid where each point can 
be labeled by 2; = {zo,zi,Z2) with Zq — zf — Z2 = —1. This applies equally well to the 
solutions we found in Section II or the Liischer-Schechter solutions discussed in the previous 
section. Any solution to (2.20) gives a p{z) and a '({^{z) from which we can generate new 
solutions p{Az) and 'ip{Az) where A is an element of 50(2,1). Note that this is possible 
because the transformation A maintains the form of the metric ds^ = dzQ — dzf — dz2, that 
is the transformation is an isometry. Thus, in general, for each solution to (2.20) there is 
actually a three-parameter family of solutions. We call this S0{2, 1) which acts directly on the 
hyperboloid, SO{2, 1)^- There is a very closely related SO{2, 1) associated with the conformal 
group in (3+1) dimensions. To understand this SO{2, 1) recall that pure Yang-Mills theory in 
(3+l)-dimensions is invariant under the 15 parameter conformal group S0{4:, 2). Working in 
a spherical ansatz, which entails picking an origin, reduces the symmetry to SO{3) x SO{2, 1). 
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The S0{3) is that of spatial rotations whose action on the gauge field A^, given by (2.2), 
is equivalent to a constant gauge transformation of the form (2.8). Since p and i/j are gauge 
invariant functions of r and t, this SO{3) leaves p and ip unchanged. We call the SO{2, 1) 
subgroup of SO{4:, 2), SO{2, 1)^. As we will see SO{2, 1)^ is composed of time translations, 
dilatations and special conformal transformations associated with time translations. We now 
discuss these two SO{2, l)'s and their relationship. 

SO(2,l)„: 

The SO{2, 1)„ is the isometry group of the hyperboloid. It takes p{z), ip{z) into p{Az), 
'ip{Az). Note that although w and r cover the hyperboloid, the coordinates r and t do not. 
From (2.21) we see that zq + Z2 = 1/r, so r and t only cover the half of the hyperboloid, 
i?"*", on which zq + Z2 > 0; see Fig. 1. For a fixed r and t we determine a z E H'^ and our 
new solution at r,t is p{h.z), 'ip{Az) even if Az is not in H~^. Thus the SO{2, 1)^ acts on the 
whole hyperboloid, even where r and t are not defined, and in this sense is not a space-time 
symmetry. 

We will decompose SO {2, 1)^ into three one-parameter transformations which together 
generate the full SO{2, 1)^. The obvious choice for the three transformations would be those 
that leave 2; = (1, 0, 0), z = (0, 1, 0) and z = (0, 0, 1) invariant. However to see the connection 
with SO {2, l)c it is more convenient to pick a different set. The first is 



At(c) = 



_1^2 _„ 1 _ 1^2 



(4.1) 



which depends on —00 < c < 00. Note that zq + Z2 is invariant under this transformation so 
— > If we use (2.21) to discover the action of the transformation on r and t we find 

r ; t^t + c , (4.2) 
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so clearly generates time translations. The second transformation is 

/ cosh a — sinh ct \ 

Ad (a) 



ii 1 
— sinh a cosh a 



(4.3) 



/ 



with — oo < q: < oo. Here zq + Z2 ^ e "'{zq + Z2) so the sign of zq + zi is invariant and 
B.^ B.^ . Acting on r and t this gives 



r — > e r ; 



t 



(4.4) 



which is the dilatation transformation. The two transformations (4.1) and (4.3) together form 
a two-parameter non-semisimple subgroup of S0{2,1)h which maps H'^ into H'^. 
The third transformation is 



/I + 1^2 -d -\d^ \ 



As{d) = 



\ 



-d 



1 



d 



-d \-\d^ J 



(4.5) 



with —00 < d < 00. For a fixed d, if (1 — dt^ — d^r^ > then the z associated with r, t is 
mapped under (4.5) into H'^ and we have 



(1 - dtf - d^r'^ ' 
t + (i(r2 - t^) 



(4.6) 



{1 - dtf - d^r"^ ' 

However if (1 — dtf — d'^r'^ < 0, then the z associated with r, t is mapped into the half 
of the hyperboloid on which zq + Z2 < 0. Note that for any d ^ 0, there are always points 
z e H'^ which are mapped into H~ . We will have more to say about the transformation 
(4.6) when we discuss SO{2,l)c- Also note that if we define the discrete transformation 
Io{zq, Zi, Z2) = {—Zq, Zi, Z2) then we see that 



As(d) = Iq K{d) lo 
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(4.7) 



In principle any solution of (2.20) which lies on the whole hyperboloid corresponds to 
a three parameter family of solutions generated by SO (2, 1)^. The solutions introduced in 
Section II depend only on zq and not on the ratio zxj zi. Thus the transformations which leave 
zq invariant do not, in this case, generate new solutions and there is just a two parameter 
family of solutions. To see this in terms of the transformations just introduced, note that 
As((i) given by (4.5) can be written as the product of a rotation about (1,0,0) times a 
dilatation (4.3) times a time-translation (4.1), 

A,((i) = i?o(/3) A.(ln(l + (i2)) A^(-_i_) . (4.8) 

Here Ro{P) is a rotation about z = (1, 0, 0) by the angle P where cos /3 = (1 — d^)/{l + d'^) 
and sin/? = 2d/{l + d'^). Since Ro{P) leaves the solutions of Section II invariant we see that 
As{d) acting on the solutions can be expressed in terms of the transformations A^, and A^. 

Except at special values of e, the Liischer-Schechter solutions are discontinuous on the 
hyperboloid along r = ±7r. Recall that r = ±7r is in H~ . The transformations At(c) and 
Aj:,{a) take — > so these transformations do not bring the cut into and therefore 
new solutions, generated by At(c) and Ao(q!), are continuous. The same can not be said of 
As{d) since by (4.8) we see that the r = ±7r line is rotated into by Rq for any (3 0. 
Thus the transformation As{d) acting on the solutions of Section III does not generate new 
solutions. (We do not call a function a solution to a differential equation if it solves the 
equation everywhere except along a line where the function is discontinuous.) The Liischer- 
Schechter solutions can be parametrized by e and tq together with the parameters a and d 
associated with dilatations and time-translations. For the values of s for which p and ip match 
at T = ±7r and there is no cut, one may think that As{d) generates new solutions. However, 
these solutions can be obtained from the original solutions by a shift of tq. 
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Because the Liischer-Schechter solutions are discontinuous along r = ±7r on the hyper- 
boloid, the transformation Rq{(3) under which r ^ r -\- (3 can not be used to generate new 
solutions since the transformation brings the cut into . However, the cut exists only be- 
cause we insist on identifying r = tt with r = — tt. If we look at the differential equations 
(2.23) we can imagine producing a solution for —oo < r < oo. In this sense we would be 
solving on a multiply covered hyperboloid and there would be no cuts. The transformation 
Rq{P) would take p(t, w), '0(t, w) to p{t + P, w), 'iIj{t + P,w). However, shifting r is equivalent 
to shifting tq (see (3.5) and (3.6)), so actually no new solutions are introduced by adopting 
this point of view. 

SO(2,l)c 

As we discussed in the beginning of this section, the full 5(9(4,2) conformal group in 
(3+1) dimensions has an SO{2,l)c subgroup which commutes with the SO{3) of rotations 
about X = 0. This group can be decomposed into three continuous transformations which 
generate the full group: 

(i) time translations: t — > t + c, a;* — > a;*. Acting on functions that depend only on r and t 
this is equivalent to (4.2). 

(ii) dilatations: t — > e~"t, —>■ e~"a;\ Again acting on functions that depend only on r and 
t on this is equivalent to (4.4). 

Thirdly, we have the special conformal transformation associated with the time translation. 
This is generated by an inversion, x^^ — > x^/x"^, followed by a time translation, followed by 
another inversion. The net transformation is 
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(iii) special conformal: 



t 



(1 - dtf - (Pr^ ' 

(4.9) 

t + d(r2 - t^) 
(1 - dtf - (Pr^ ' 

Note that along the light cone, (1 — dtY' — d'^r'^ = 0, the denominators in (4.9) vanish and the 
transformed coordinates blow up. As you cross this light cone and (1 — dtY — d'^r'^ changes 
sign, the transformed coordinates change discontinuously from say — oo to +oo. This can 
induce discontinuities in the transformed solution unless the original solution is particularly 
well behaved as a;* — > ±oo and t — > ±oo. 

Acting on functions which depend only on r and t, the special conformal transformation 
given above is 



{1 - dtf - d'^r^l ' 

t + d{r^ -e) ^ ' 



(1 - dty - d^r"^ ' 

which is the same as (4.6) only when (1 — dt)'^ — d'^r'^ > 0. The S0{2, l)c is a space-time 
symmetry. It takes fields defined on x, t and generates new fields on x, t using the old values 
of the fields. This is distinct from SO{2, which can take values of the fields on H~ and 
bring them into H~^. For space-time functions p{r,t), '^{r,t) we can still understand how 
SO{2, l)c acts by looking at the hyperboloid. For each r, t there is a 2; e , so given p(r, t) 
and '^{r,t) we know p{z) and '^{z) on H'^ . Now define new fields p{z) and '^{z) by 

p{z) = {P^'^\ (4.11) 



and similarly 



V'(^)=^f!'^\ (4-12) 
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Thus p and t/^ are defined everywhere on the hyperboloid. We now define p'{z) = p{Az) and 
iIj'{z) = iIj{Az) with A e SO {2, from which we can infer p'(r, t) and '(/''(r, t). For A = Aj.{c) 
of (4.1) it is immediate that p'(r, t) = p(r, t + c) and ip'{r, t) = 'ip{r,t + c); for A = Ao{a) of 
(4.3) it is immediate that p'{r, t) — p{e~°'r, e~°'t) and il^'ir^ t) — 'i(j{e~'^r, e~"t). For A = As{d) 
of (4.5), for which Az can be in H~ , we use the fact that z — > —z is equivalent to r — > — r in 
(2.21) to show that p'{r,t) — p{r',t') and il:'{r^t) — '0(r',t') with r' and t' given by (4.10). 

In short the action of S'0(2, 1)^^ and S0{2,l)c can be summarized as foUows. For a 
solution which is defined on the whole hyperboloid, SO{2, 1)^ defines a new solution p{Az), 
'ip{Az) with A e SO{2, 1)„. The SO{2, l)c is a space-time symmetry and only refers to p{z) 
and ip{z) defined on 2; e H'^ . To find the action of SO{2,l)c you discard p{z) and ip{z) 
for z & H~ and replace them with p{—z) and '4>{—z) giving p and defined on the whole 
hyperboloid (see (4.11) and (4.12)). The SO{2, l)c gives new solutions p'{z) = p{Az), il)'{z) = 
■ip{Az) with A e S0{2, 1)^. For time translations and dilatations the two S0{2, l)'s act in 
precisely the same manner. However, for A's given by (4.5) they are different transformations. 
For a fixed d and r,t such that (1 — dtY — d^r^ > 0, the two transformations agree. However 
for (1 - dt)'^ - d'^r^ < they disagree. 

In using (4.11) and (4.12) we should note that the boundary of and H~ is the intersec- 
tion of the plane zq + Z2 = Q with the hyperboloid. For p{z) and '^{z) to be continuous on the 
hyperboloid we require p{zq, zi, -zq) = p{-zo, -zi,zo) and ip{zo, zi, -zq) = ip{-zo, -zi,zo). 
If these conditions are not met, then p and ip are discontinuous along zo + Z2 = 0. In this case 
p{Asz) and '(/'(A^^) viewed as functions of r,t are discontinuous and the special conformal 
transformation does not produce new solutions. The solutions of Section II are functions of 
Zq only, but these functions are not even, so the special conformal transformation acting on 
these solutions produces ripped functions. Similarly no new solutions are produced by acting 
with special conformal transformations on the Liischer-Schechter solutions. 
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V. TOPOLOGICAL PROPERTIES OF CLASSICAL SOLUTIONS 

In Section II we found a new class of Minkowski space solutions to the Yang-Mills equations 
and in Section III we described another class of solutions discovered by Liischer and Schechter. 
Both of these classes have the following general properties. At a large negative time, —Tq, the 
energy density of the solution is localized in a thin spherical shell with radius of order Tq, and 
the shell is collapsing. Imagine a 2-sphere of radius R, which we call Sr, which is concentric 
with the energy shell. For Tq » R the energy density of the classical solution is initially 
localized far outside of Sn so the region inside this sphere is very close to pure gauge. At a 
time t ~ —R the energy front reaches Sr. Energy flows in, bounces back at time t ^ 0, and 
then leaves the region at time t ~ i?. At large positive time, Tq, the fields inside Sr are once 
again pure gauge. But this final vacuum configuration need not coincide with the initial one. 

The Local Winding Number Change 

In the Ao = gauge, the initial and final vacuum configurations inside Sr are given by 
pure gauges of the form (2.7). Such configurations inside the sphere can be characterized by 
the quantity 



where the integration is over the interior of Sr. We will refer to this expression as the local 
winding number. A pure gauge configuration in the interior of Sr with U — 1 (or any constant 
element of SU{2)) on the boundary defines a map from compactified three space - the three 
sphere - into the gauge group. In this case ^'(i?) is an integer, the winding number, which 
characterizes the homotopy class of the map. If U is not required to equal unity (or a constant 




(5.1) 
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element of SU{2)) on Sr, ^{R) will not in general be an integer. Furthermore if we do not 
restrict U on Sr, then ^{R) will not be additive on the products of two Vs. 

The solutions which we consider excite the interior of Sr only for —R ^ t ^ R. Before 
and after this period the field inside the sphere is pure gauge. We are interested in calculating 
the winding number change of these pure gauge configurations. Working in the = gauge 
we are free to make time-independent gauge transformations. We use this freedom to set 
U = 1 at points inside of Sr at very early times. Therefore ^'(-R) = at these early times. 
Inside of Sr there is no gauge freedom left and the form of U inside of Sr at late times is 
now determined by the classical solutions. The quantity we are after is (5.1) evaluated with 
U determined by the late time behavior of the solutions. 

The prescription we have just given for determining the local winding number change 
associated with a solution is in fact gauge invariant. Consider a solution A^(x, t) which has 
the property that at very early and very late times, is pure gauge for |x| < R. Go into the 
^0 = gauge. At very early times Aj = ^UEdjUl and at very late times Aj = ^U^djUl for 
|x| < R. Inside of Sr make a gauge transformation by C/1 so that in the far past Aj = 0. Now 
at late times Aj = ^UdjU^ with U = UlUi^. We then evaluate z^(-R) given by (5.1) for this 
U. If we had started with a gauge transform of A^, say and then passed to the = 
gauge we would have discovered that Ue goes into VUe and goes into VU^ where V is 
time independent. Thus U — uIUl is gauge invariant. We are evaluating the local winding 
number of the difference between Ue and which is uIUl. Note that the local winding 
number of C/lC/i is not the difference between the local winding numbers of C/^, and Ue since 
z^(i?) is not in general additive. 
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We now return to the spherical ansatz where U has the form (2.8). Substituting this 
form into (5.1) gives the following expression for the local winding number inside of Sji in 
terms of the gauge function Vl{r) which characterizes the late time pure gauge inside of Sr: 



zy(i?) = J- /"^dr [l-cosn(r)] . 

27r Jo dr 



(5.2) 



Note that the expression for the local winding number (5.2) which arises from a (3+1)- 
dimensional discussion is different from the expression for the local winding number, uj, which 
one might naively infer from the (l+l)-dimensional U (1) gauge theory vacuum given by (2.9), 



1 



R 



dfi(r) 



(5.3) 



The right hand sides of (5.2) and (5.3) coincide, however, if f2(r) is an integer multiple of 27r 
at r = and r = R. This illustrates that the definition (5.1) of the local winding number is 
somewhat arbitrary. Indeed, you might pick i'{R) to be any function of U which coincides 
with (5.1) when U is restricted to a constant on Sn. 

We now express the local winding number (5.2) of the final vacuum inside of Sr as t — > oo 
in terms of finite-energy classical solutions. Using (2.9) and (2.11) we see that (p{r,t) — > Q{r) 
as t — > +00 for any fixed r. Thus, (5.2) gives 



u{R) = ^ [ <f{r, oo) - sin (p{r, oo) ] 



--R 



r=0 



(5.4) 



Now the finite energy condition (2.19a) gives p{0,t) = 1 for all t. In fact, it can be shown 
that all finite energy solutions obey drip{r, ^)|^^q = 0, from which (2.15) implies dt(p{0, t) = 0. 
Thus we have that at r = 0, <^ stays locked down at its initial value, which we have taken 
to be zero, so ip{0, t) — 0. Recall from our discussion following (2.9) that we must require 
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0,{r = 0) to be equal to zero for pure gauges in (3+l)-dimensions. We now see that this 
condition is dynamically imposed on all finite energy solutions. Therefore equation (5.4) can 
now be written as 

^'(i?) = r <^(i?, oo) — sin(^(i?, oo) 1 . (5.5) 

Note that ^'(-R) is intimately related to the change in the phase of the field x at produced 
by the passing energy front. 

We are interested in evaluating (p{R, oo) in the ao = gauge. In the late-time limit, 
t — > oo, the fields are pure gauge for any finite r, so dr(fi{r, oo) = ai(r, oo). Thus 

(p{R, oo) = / dr ai{r,(X)) 
Jo 

dt / dr doai{r,t) (5.6) 

-oo Jo 

-I POO pR 

= --j dt dre^^f^, , 
where the last expression is manifestly gauge invariant. Using (2.13) we have 

V?(i?,oo) = -2 r dt r dr , (5.7) 

J-oo Jo f 

which can be expressed in terms of the hyperboloid variables w and r as 

ifiR, oo) = -2 / dwdr ti^!!lll , (5.8) 
Jhj^ cos^ w 

where Hn is the part of for which r < R. From (2.21) we see that the fixed R contour is 
given by the intersection of the plane zo + Z2 = 1/R with the hyperboloid. In terms of w and 
T, from (2.22), we have that at fixed R, w and r obey 

cos T = — cos + sin . (5-9) 
R 
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The region Hr is shown in Fig. 7. We can now write (5.8) as 

/■«^max rr+(.w,R) , / N 

cp{R,oo) = -2 dw / dr y' , (5.10) 

J-n/2 J-T+{w,R) cos W 

where T+(t(;, R) is the positive value of r which solves (5.9) for fixed R and tu, and tUmax obeys 

i?^ — 1 

sin Wmax = ^2 _^ 2 ' ''^""^"^^ 

It is now straightforward to evaluate (5.10) for any solution in the spherical ansatz. 
The Local Winding Number Change Produced by the Solutions of Section II 

Here we discuss ^{R) given by (5.5) for the solutions of Section II. There is no reason 
to expect z^(-R) to have an integer value at any R and in fact i'{R) diverges for the solutions 
of Section II as i? — > oo. In evaluating (5.5) we use (5.10) and note that in this case ij) is 
independent of r, so 

/•"'max J^l \ 

(p(R,oo) = -4 dw T+iw,R) . (5.12) 

C0S2 W 

For w near — 7r/2, the integrand is integrable because ipi^w) goes to zero fast enough; see 
(2.28b). For finite i?, Wmax is less than 7r/2 and accordingly ip{R,oo) and i^(-R) are finite. 
In Fig. 8 we show z^(-R) versus R for a typical solution. Note that i^iR) does not appear to 
have a limit as R goes to infinity. This is indeed the case which can be seen from (5.12). As 
i? ^ oo, Wmax — ^ 7r/2 and t+{w, R) ^ n/2 — w. Since ip{'n/2) ^ the integral diverges at its 
upper limit. 

The Local Winding Number Change Produced by the Liischer-Schechter Solutions 

The local winding number ^'(-R) of a Liischer-Schechter solution approaches a fixed value 
as i? — > oo. To see this we again evaluate (5.5) using (5.10) for fixed R. However in this case 
from (3.2b), 'ip{w,T) = ^q{T)cos'^w so 



max 



<fi{R,oo)= dw {q[-T+{w,R)]-q[T+{w,R)]} . (5.13) 
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Recall from the discussion in Section III that q{T) can be viewed as the coordinate of a particle 
moving in a potential U{q) with "energy" e as expressed by (3.4). Besides £, the solutions 
are characterized by tq, the "time" when the particle is at a turning point of U{q). If tq = 0, 
then q{T) is an even function and (5.13) vanishes for all R. To have non-trivial topological 
properties we must look at tq ^ 0. In Fig. 9 we give examples of versus R for tq = 1 
and e = 0.7 and 2.0. Note that as i? — > oo, v{R) appears to approach a limit. To obtain a 
simple expression for this limiting value, note as we mentioned just before, that as i? — > oo, 
'ii'max 7r/2 and Tj^{w^ R) — > 7r/2 — w. Thus as R goes to infinity, (5.13) gives 



In Figs. 10 and 11 we show the asymptotic winding number v = limi^^oo i^{R) which is 
obtained from (5.14). We plot both u versus e for fixed tq and v versus tq for fixed e. In the 
latter case the asymptotic winding number is periodic in tq. This period is exactly the period 
of a particle moving in the potential U (q) with energy s. 

The Local Topological Charge 

We first discuss the local topological charge of any solution in the spherical ansatz. We 
consider the local topological charge defined as 



where the spatial integration is over the interior of the sphere of radius R. This quantity is 
manifestly gauge invariant. The topological charge, Q, given in (1.1) is the limit of (5.15) as 
R ^ oo. Since the topological charge density is the divergence of a current, the topological 




(5.14) 




(5.15) 
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charge can be expressed as a surface integral. In the spherical ansatz we can write the following 
expression^^'-'^'^ for Q{R), 

dt / drd^f , (5.16) 

-oo ^0 



where the two dimensional current j^^ is 



27r 



ay - Re + 2^ ( X*DyX - xiD^x)' 



(5.17) 



Integrating (5.16) gives 



Q{R) = / dr f 
Jo 



t=oo 



t=-oo J-oo 



+ / dt f 



r=R 



r=0 



(5.18) 



The first term we recognize as the change in winding number J>'(i?) given by (5.5) and the 
second term is the net flux through R. 

For solutions to the equations of motion ^^j'^ takes a simple form. In terms of the gauge 
invariant variables p and i/j we have 



Q{R) 



2tt 



-R 

dt I dr 



'-00 ^0 



(5.19) 



which upon using (5.7) gives 



Q{R) = —ip{R, oo) + J^R 



(5.20) 



where the first term in this decomposition of Q{R) differs from i'{R) by ^ sin[<^(i2, oo)] and 
the second term, JFr, is defined as follows: 



-1 /-oo nR 

^R=^J dt dr d^ir^dyiP) 



(5.21) 
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Here again ry'*^ = diag(-l, +1). Since ds^ = {dr"^ - dt^)/r'^ = {dw"^ - dr"^)/ cos^tu it follows 
that 

^R=^f dwdrdair^d^iP) , (5.22) 
where a, P = t,w. We can express Tji as an integral along the boundary of Hji: 



R 



1 r 



dr 



dw dr dr 



(5.23) 



W=Wb{T,R) 



where Wb{T, R) is the value of w determined by (5.9) for fixed R and r. In deriving (5.23) we 
have used the fact that for finite energy solutions §^ I ,„ = 0, so there is no contribution 

ow \w=—'k/2 ' 

along the lower w = — 7r/2 boundary of Hn as seen in Fig. 7. We will shortly evaluate J^r 
for the two classes of solutions we have discussed and obtain Q{R) from (5.20). 

We can also obtain a simple expression for lim/^^oo ^R- Note that as i? — > cxo, — > 
7r/2 + T for T < and tu^ — > 7r/2 — r for r > 0. In terms of the variables w± = w ±t, (5.23) 
implies 



lim Tji= [ 
-R^oo 27r J_ 



dw^ 

3n/2 



dw. 

/2 



dw_ 



ui-|_=7r/2 



(5.24) 



1 



TT 



ip{w = 7r/2,T = 0) 



where we used the fact that for all finite energy solutions il){w = —'k/2,t) = 0. For so- 
lutions with the property that (fi{R,oo) has a finite limit as i? ^ oo, (5.10) implies that 
'ip{w = 7t/2,t = 0) = 0. In this situation we obtain 



<5 = 77- Jim (fi{R, oo) . 

ZTT it— »oo 



(5.25) 
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The Local Topological Charge of the Solutions of Section II 

For the solutions of Section II where ij) depends only on tu, we have from (5.23) that 



For a fixed value of R we can evaluate this expression given a numerically generated solution 
il){w). To get Q{R) we use (5.20) with Tr given by (5.26) and <^(-R, oo) given by (5.12). An 
example of Q{R) as a function of R is shown in Fig. 8. Again Q{R) has no limit as i? — > oo 
since Lp{R,oo) does not approach a limit. 

The Local Topological Charge of the Liischer-Schechter Solutions 
Here again ip{w,T) = ^^(t) cos^ Thus from (5.23) 



Now if To = 0, g(T) is an even function and J^r — 0. Thus for tq = 0, Q{R) — 0. For tq ^ 0, 
in general !Fr will not vanish. We can evaluate !Fr for a given solution and then with the use 
of (5.20) we can determine Q{R). Examples are given in Fig. 9. 

For the Liischer-Schechter solutions, <^(-R, oo) has a limit as i? — > oo as can be seen from 
(5.14). Thus (5.25) applies and we can obtain Q which we see is not an integer. Figures 10 
and 11 show the asymptotic topological charge for typical Liischer-Schechter solutions. 
Integer Topological Charge? 

You may think that Q = limi^^oo Q{R) should be an integer because of general topo- 
logical arguments. However, this is not so, essentially because finite energy solutions are to 
be found at arbitrarily large radii for arbitrarily large times. The argument which leads to 
integer values of Q is as follows. Consider a region of space-time which contains non-zero 




(5.26) 




(5.27) 
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energy and imagine that outside this region the energy density is zero. You can surround 
the region of space-time by a three dimensional surface which is topologically S^. On this 
surface the gauge field is pure gauge, i.e. — i/g Udf^W . Thus we have a map from 
into the gauge group which is characterized by an integer. The topological charge, integrated 
over the space-time region in question, gives this integer. For non-zero energy solutions to 
the equations of motion, it is not possible to surround the energy density by a three-sphere in 
space-time because the energy density is moving out to infinity at early and late times. We 
do not expect and do not find integer values of Q. 

VI. DISCUSSION 

Our ultimate aim is to relate the classical solutions discussed in this paper to physical 
processes. In a companion paper ^ we analyze fermion production in the background of a 
field which locally changes winding number. We work with a (1-1-1) -dimensional analogue. 
We show that non-integer change in winding number is associated with quantum mechanical 
fermion number production where the change in winding is the expectation of the number of 
fermions produced. Therefore we believe that if we could produce a field configuration like 
one discussed in this paper, it would be associated with fermion number violation. Here we 
are ignoring the back-reaction that the fermions have on gauge fields. 

We are then led to seek the relationship between the classical solutions and the gauge 
boson quanta of the real world. If the classical solution represents a coherent quantum state, 
then we need to understand the overlap of this state with the few particle quantum state of an 
accelerator beam if we are going to use our methods to estimate rates for fermion violation in 
actual experiments. We plan to pursue this as well as the question of the relative weight these 
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coherent states have in a high temperature density matrix which should tell us the relevance 
of these solutions to high temperature processes. 

We are intrigued by the fact that these solutions change the winding number of local 
regions of space by non-integer amounts and create non-integer topological charge. Using 
a solution to the classical field equations you can excite an arbitrary large region of space 
and discover that after the energy has dissipated from the region, the winding number of 
the region will have changed by a non-integer amount. The implications of this for the full 
quantum theory, we have yet to discover. 
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APPENDIX 

In this appendix we investigate the analytic structure of solutions (2.24) near the end 
points tUinit = — 7r/2 and w^n = +7r/2. From (2.21) and (2.22) we can see that w-mit corre- 
sponds to r = for fixed time t, we_n corresponds to r = oo for fixed time t, and that t = ±oo 
for fixed r also corresponds to w-mit- Recall from (2.19) that for finite energy solutions p is 
near one and ijj is near zero for w near Winit- Writing w = winit + e, from the differential 
equations (2.24) we can infer that for e << 1, 

p = l + 0(e2), (Ala) 
iP = 0{e^) . {A.lb) 

For fixed t and small r, e ~ r so (A.l) gives 

p = 1 + 0{r^) , {A.2a) 
^ = 0{r^) , {A.2b) 

while for fixed r and large \t\, e ~ so (A.l) gives 

p = 1 + 0{t-^) , (A.Sa) 
V; = 0{t-^) . (AM) 

Note that (A. 2) is more stringent than just the finite energy bounds (2.19). 

Working near w^n when r is large, the finite energy condition places only weak restrictions 
on p and ip. However, to solve (2.24) near wnn = 7r/2 we must have p ±1,0. Writing 
w — Wfln — e, the asymptotic form for small e in the case that p — ±1 becomes, as in (A.l), 
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p = ±1 + 0{e^) , (AAa) 
^ = 0{e^) , {AAh) 



whereas if p ^ we get 



p = a e"^ , {A.5a) 
^ = ^^ + be^ , (AM) 

where rj and ^ are as yet undetermined. We examine the p — > case. Substituting (A. 5) in 
(2.24b) gives ^ = 2r]. After substituting (A. 5) into (2.24a) we find two complex solutions for 
r] when b ^ 0, 

V± = l l±iy^3 + 16V^ . (A6) 

Complex solutions imply oscillatory behavior; however, since the equations are non-linear (for 
small p the 1/p^ in (2.24a) is important) we cannot superimpose these solutions to form real 
ones. It turns out that we can linearize (2.24a) in p in the following manner. To leading order 
in e, (2.24b) implies that 

^ - . (A7) 

p^(e) e 

where a prime denotes differentiation with respect to e. Using this form for ip' in (2.24a) we 
have to leading order in e 

,"+i±M, = o. (AS) 

The general solution to (A. 8) is p = a^e'^+ + a-e^- , where a± are arbitrary constants and 'q± 
are given by (A. 6). We can now form a real solution by appropriate linear combinations. We 
write 

p = cos[c Ine + B] , {A.9) 
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where c = •\/3~+T6V^/2 and A, B are arbitrary constants. Substituting (A. 9) back into (A. 7) 
and integrating gives 



^ ^ cos(2c Ine + 2B) ^ 2c sin(2c Ine + 2B) 



(AlO) 



1 + 4c2 1 + 4c2 

The solutions found numerically in Section II appear to possess this behavior. Only two of 
the constants A, B and ipo in (A. 10) are independent due to the two parameter nature of 
finite energy solutions. 

We now prove the assertion that for finite energy solutions of (2.24), if ip is non-zero 
then it is monotonic, that is to say ip is non-increasing or non-decreasing. We prove this 
by contradiction. From the finite energy boundary condition (A. lb), iIj{w) must vanish at 
the point Winit = —n/2. Let w* be the first point greater than Wjnit for which ip ceases to 
be monotonic. This means ip is either non- increasing or non-decreasing from tUinit to w*. 
Assume the latter (which is consistent with Fig. 4); the former case may be handled in a 
similar manner. Then w* is a relative maximum of ip{w), i.e. '^'{w*) = and ip{w*) > ip{w) 
for points w slightly to the left and right of w*. Since ip is non-zero with vanishing derivative 
at w*, (2.24) implies that p{w*) = 0. We write w = w* + e and expand p and ip in powers of 



p = aie + a2e^ + a^e^ + . . . , (^.lla) 
^ = bo + b2e^ + b3e^ + b4e^ + ... . (A.Ub) 

For ip to have a relative maximum at w*, the leading non-trivial e dependence must be an 
even power of e with a negative coefficient. We will now show that this can not happen. 

Assume that ai ^ 0. We will relax this assumption in a moment. Note that bo > since 
ip started out zero and never decreases until w*. Using (A. 11) in (2.24b) we find 
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b2 = 







(Al2a) 



80263 
2ai 



+ 



(yl.l26) 



2 cos^ w* 



while 63 is left undetermined. We are not interested in the case in which 63 ^ 0, since the 
leading e behavior would then be an odd power in e. So if 63 = 0, then (A. 12b) gives 64 > 
since bo > 0. This is a contradiction since the leading epsilon dependence must be an even 
power of e with negative coefficient. If ai had vanished, let a„ be the first non-zero coefficient 
of (A. 12a). Then similar reasoning would give ip = bo + |62n+2|e^"'"'"^ + . . ., which again is a 
contradiction. Therefore, iJj{w) is monotonic over the whole w-range. 

A corollary of this result is that the only finite energy solutions that interpolate between 
p = +1 and p = — 1 as varies from tUinit to w^n have i/j equal to zero. This can be seen as 
follows. Finite energy solutions start at p = 1 and •0 = 0. From (A. 4), if p does not asymptote 
to zero as — > wa-n then i(; must return to zero. This can not happen for non-vanishing i/j 
which never decrease. 
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FIGURE CAPTIONS 

Fig. 1 The hyperboloid ZQ—zf — Z2 = —1. The r,t coordinates cover the half of the hyperboloid, 
H~^, above the zo + = plane. We also show a typical fixed-r contour Cr. 

Fig. 2: The function p{r]) for the initial conditions p(0) = 0.2 and p(0) = —0.881146. For a given 
p(0), the velocity p(0) must be fine-tuned to ensure the finite energy condition p — > 1 as 
?7 — > — oo. 

Fig. 3: A typical p{w) for the solutions of Section II. We have chosen the initial conditions 
p{wi) = 0.99995, p'{wi) = 0.01, ip{wi) = -0.0001 and ip'{wi) = 0.02. To avoid the 
singularity in (2.24) at —n/2 we have taken Wi = —n/2 + 10~^. In all graphs related to 
Section II, we will use the initial conditions stated here. 

Fig. 4: A typical ■ijj{w) for the solutions of Section II. 

Fig. 5: Profiles of the energy density, e{r,t), times in units of Syr/gr^ for a sequence of times 
for a typical solution of Section II. 

Fig. 6. Profiles of the energy density, e(r, t), times in units of Sir/g'^ for for a sequence times 
for a Luscher-Schechter solution with £ = 0.2. 

Fig. 7. The regions H'^ and Hu in the w-r plane. 

Fig. 8: The local winding number i'{R) and the local topological charge Q{R) vs. R for a typical 
solution of Section II. 

Fig. 9: The local winding number ^{R) and the local topological charge Q{R) vs. R for the 
Luscher-Schechter solutions with two values of e and tq — 1. 
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Fig. 10: The asymptotic winding number v and the topological charge, Q, for the Liischer- 
Schechter solutions as a function of e with tq = 1. 

Fig. 11: The asymptotic winding number v and the topological charge, Q, for the Liischer- 
Schechter solutions as a function of tq with e = 0.7. 
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